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Luttinger theorem for the strongly correlated Fermi liquid of composite fermions 
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While an ordinary Fermi sea is perturbatively robust to interactions, the paradigmatic composite- 
fermion (CF) Fermi sea arises as a non-perturbative consequence of emergent gauge fields in a 
system where there was no Fermi sea to begin with. A mean-held picture suggests two Fermi 
seas, of composite fermions made from electrons or holes in the lowest Landau level, which occupy 
different areas away from half Filing and thus appear to represent distinct states. We show that in 
the microscopic theory of composite fermions, which satishes particle-hole symmetry in the lowest 
Landau level to an excellent degree, the Fermi wave vectors at Filing factors h' and 1 — ly are the 
same, and are generally consistent with the experimental Fndings of Kamburov et al. [Phys. Rev. 

Lett. 113 , 196801 (2014)]. Our calculations suggest that the area of the CF Fermi sea may slightly 
violate the Luttinger area rule. 

PACS numbers: 73.43.-f 


A fundamental property of a Landau Fermi liquid is 
captured by Luttinger’s theoreuP, according to which the 
volume occupied by the Fermi sea, appropriately deFned, 
remains invariant when the interaction is switched on, so 
long as no phase boundary is crossed. A violation of this 
theorem signiFes non-Fermi liquid behavior, which has 
motivated invest igation^^H^ of its applicability for various 
strongly correlated systems, such as high temperature 
superconductors. This article investigates the Luttinger 
theorem for an exotic emergent Fermi sea. 

When two-dimensional electrons are subjected to a 
strong magnetic Feld, they exhibit the phenomenon of 
the fractional quantum Hall effect (FQHE)P, which is un¬ 
derstood in terms of topological particles called compos¬ 
ite fermionJ^inni. Halperin, Lee and Reacff^and Kalmeyer 
and Zhan^I31 theoretically predicted that at Landau level 
(LL) Filing factor u = 1/2, the external magnetic Feld 
is canceled, in a mean Feld (MF) approximation, by the 
emergent gauge Feld carried by composite fermions, and 
they form a Fermi sea. The composite fermion (CF) 
Fermi sea is special in the following sense. Ordinarily, 
we begin with the perfect Fermi sea of non-interacting 
fermions and then ask how interactions degrade or de¬ 
stroy it. In contrast, interactions are fully responsible 
for creating the CF Fermi sea (CFFS) in a system of 
electrons conFned to the lowest LL (LLL) where, orig¬ 
inally, there was no Fermi sea, and, in fact, no kinetic 
energy and no composite fermions. The very existence 
of the CFFS thus is a manifestation of strong correla¬ 
tions. The essential validity of the CFFS has been con- 
Frmed in extensive detail in many experiment^^^^^, and 
it also dovetails with the prominently observed sequences 
of fractions at v = nj (2n ± 1^. 

Kamburov et have recently made accurate mea- 
surements of the CF Fermi wave vector through commen- 
surability effects in the presence of a periodic modulation. 
They have observed more commensurability oscillations 
than before, and thus provided the most detailed conFr- 


mation of the CFFS state to date. Furthermore, the un¬ 
precedented accuracy of their measurement has revealed 
an intriguing puzzle. For electrons conFned to the LLL, 
one can take two exactly equivalent starting points: One 
can deFne the problem in terms of either electrons at v 
or holes at 1 — v. One can then go ahead and composite 
fermionize either electrons or holes to produce what we 
will label ®CFs or ^CFs, which experience an effective 
magnetic Feld given by H — 2^00, where = hc/e 

and p is the density of composite fermions. (All CF quan¬ 
tities are marked by an asterisk.) For v ^ 1/2 the ®CFs 
and ^CFs have different densities, producing, in the MF 
description, different Fermi wave vectors for the fully spin 
polarized CFFS state: 

MF for ^CFFS : = ^4^ k^i = 

MF for ^CFFS : k^ = ^4^ k^i = v/2(l - i^) 

where i = ^^hc/eB is the magnetic length, and the elec¬ 
tron and hole densities are given by pe = pu = n/{27rP) 
and ph = pi-u = (1 —z^)/(27r^^). The CFFS thus appears 
to have a split personality. This raises many interesting 
conceptual questions. At a given n, do the ®CFFS and 
^CFFSs represent two distinct states, or are they dual 
descriptions of the same state? If the former is true, 
then which of these two, if either, occurs in real sys¬ 
tems? If the latter is true, how does one reconcile the 
seemingly incompatible consequences of the MF picture, 
and how does one understand the violation of the Lut¬ 
tinger theorem for at least one of the two descriptions? In 
either case, what is the role of particle-hole (p-h) sym¬ 
metry in the LLL? Finally, how do we understand the 
remarkable Fnding of Kamburov et al^ that the mea¬ 
sured value of the CF Fermi wave vector is consistent 
with that expected from the smaller Fermi sea, namely 
kpi = min[V^, y^2(I — u)]? 

These observations have motivated two striking the¬ 
oretical proposals that lead to experimentally testable 
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predictions. Son has proposed^ that viewing the com¬ 
posite fermion as a Dirac fermion allows a p-h symmet¬ 
ric description of the FQHE and the CFFS. Barkeshli, 
Mulligen and FisheJ^ have taken the experimental ob¬ 
servations to imply that the ®CFFS and the ^CFFS are 
distinct states of matter and a spontaneous breakdown 
of p-h symmetry within the LLL selects one of them. 
Within the C hern-Simons (CS) formulation of composite 
fermion^ISmi, the MF Fermi wave vector is not ex¬ 

pected to change to all orders in a perturbative treatment 
of the Coulomb and the gauge interactions, suggesting 
that the ®CFFS and ^CFFS are perturbatively discon¬ 
nected, i.e. are topologically distinct, for any 7 ^ 1 / 2 , 
and, by extension, also for u = 1 / 2 . The CS formula¬ 
tion, however, does not implement the LLL constraint , 
and hence does not satisfy the p-h symm etry, as has been 
stressed elsewhere in the literatur d^^ l ^^ i 

We determine CFFS area using a different theoretical 
formulation of the CF paradigm, na mely the m icroscopic 
wave functions of composite fermion^^EHHM] xhis the¬ 
ory (i) is explicitly restricted to the LLL; (ii) satisfies p-h 
symmetry; and (iii) does not assume, a priori, any spe¬ 
cific value for k^£. We show that k^£ defined from Friedel 
oscillations in the pair-correlation function has the same 
value for states at u and 1 — related by p-h symme¬ 
try. We explicitly calculate kp£ for filling factors in the 
vicinity of z/ = 1 / 2 . 

We define the Fermi wave vector through the Friedel 
oscillations in the pair correlation function, for which we 
take the fornP^ 

g{r) = 1 +sm{2kpr + 9) (1) 

where ^4, a, k^ and 0 are fitting parameters. We denote 
the particle coordinates by either rj or Zj = Xj — 
and set = 1. This form is motivated by the observation 
that for noninteracting fully-spin polarized fermions in 
two dimensions at 5 = 0 , the oscillatory part of g{r) for 
large rk^ is given by ( 4 / 7 rr^/cp) sin(2/cFr). Let us consider 
a wave function (j)j, for a uniform density state at filling 
factor V. Its pair correlation function is given by 

9v{r, r') = p-‘^{4)^\¥{r)¥{r')^{r')^{r)\4>^) ( 2 ) 


where |1) is the state with the LLL fully occupied. Sub¬ 
stituting into the expression for the pair correlation func- 
tion and noting (l|/(c„,c 1 „)|l) = ( 0 |/(c„ -)> ^ 

Cm)| 0 ) produces the relation 

5 i_^(r, r') = p'{‘^^{4>^\^{r)^{r')¥{r')¥{r)\(i)^) (5) 
In terms of the LLL projected delta functioiP212il 


5{r,r') = ^exp 



r'P - 2/* 


+ 


(6) 


which satisfies 6{r,r') = [(5(r',r)]*, we have 

u5{r',r), and = {1 — u)5{r,r'). 

Straightforward algebra gives the relation (assuming 
thermodynamic limit and translational invariance, and 
setting r' = 0 ) 


where we have assumed the same magnetic lengths for u 
and 1—u. For r ^ 1, this reduces to = l + (z//(1 — 

u))^{gjy{r) — l). The important point is that an oscillatory 
term sm{2k^r) in gjy implies identical oscillatory behavior 
for gi-u^ indicating that the states at u and 1 — z^ have 
the same k^£. The “exact” kp£ is thus the same at u 
and 1 — z/, and is independent of whether the problem is 
formulated in terms of electrons or holes. 

We next determine the value of k^£ in a micro¬ 
scopic calculation from the oscillations in g{r) follow¬ 
ing Refs. I29l^ Being an equal time correlation func¬ 
tion, g{r) can be evaluated from the knowledge of the 
microscopic wave functions for the ground state in the 
vicinity of v = 1 / 2 . We concentrate on the fractions 
v = nl{2n ± 1 ) which approach v = 1/2 in the limit of 
sufficiently large n. The microscopic Jain wave functions 
for these states are given b}El 


N 

^n/(2n±l) = ^LLL n (zj 2 ^/ 0 ) (8) 

j<k=l 


1 r ^ ^ . 

= Vl / n • • • Tat) JJ 4'^(rfc)|0) (3) 

■ i=i k=i 

where 0 zy(ri, • • • r^v) is the real space wave function, 

electron annihilation op¬ 
erator in the LLL and T’l’(r) is the corresponding elec¬ 
tron creation operator, with the single particle LLL wave 
function defined as grn = {27r2^m\)~^^‘^z^ ex.p[—\z\‘^/4]. 
We can similarly define the pair correlation function for 
electrons at 1 — z/, with 

I r ^ ^ . 

/ n • • Tiv) tJ ^(rfc)|l) (4) 

‘ j=l k=l 


where 5 lll denotes LLL projection and is the wave 
function of n filled LLs, with = [4>^]*. At first it may 
appear that the above mentioned dichotomy is present 
also in the microscopic theory of composite fermion, as 
we illustrate by considering the fraction z^ = (Sim¬ 

ilar considerations apply to z^ = n/( 2 n + 1).) According 
to the CF theory, there are two ways of constructing a 
FQHE state at this fraction: (i) As the electron partner of 
the z/* = n integer QHE (IQHE) of ^CFs in positive 5*, 
with wave function given by Cp-h^LLL Yljcki^j -Zk)‘^^n 
where Cp_h represents p-h transformation, (ii) as the 
z/* = n + 1 IQHE of ®CFs in a negative 5*, with wave 
function given by 'PLLhYijKki^j “ Zk)^[^n+i]*■ A pri- 
ori, these appear to represent two distinct FQHE states, 
and one may ask which one applies to the real system. 
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However, explicit evaluationhave demonstrated the 
nontrivial result that these two descriptions represent the 
same state. They predict identical quantum numbers for 
the ground state and the excitations (see Supplemental 
Materials and there is an almost prefect over¬ 

lap between the two wave functions wherever it has been 
evaluated (e.g. for the 10 particle 2/3 state, the two 
wave functions have overlaps of 0.996 and 0.994 with the 
exact Coulomb stat^^. We have evaluated the ^(r)’s 
of 4 ^(n+i)/( 2 n+i) and T^/(2n+i) and found them to be 
related by p-h symmetry to a very high accuracy (See 
S^PI). The wave functions T^/( 2 n±i) produce, in the 
limit of n ^ oo, the same CFFS from either side, because 
4>(H* = 0) is reaP^. Further, Rezayi and Haldan^^have 
directly constructed the wave function for the CFFS on a 
torus and found that, for N = 16 particles, it has has an 
overlap of 0.9994 with its hole conjugate, and 0.9925 with 
the exact p-h symmetric Coulomb ground stat^^. The 
degree to which the microscopic wave functions of the 
CF theory satisfy the p-h symmetry may seem surpris¬ 
ing, but is a byproduct of the fact that these are excellent 
approximations of the exact Coulomb states in the LLL 
which satisfy p-h symmetry exactly. 

The understanding of FQHE at = (n-hl)/(2n + l) as 
z/* = n + 1 IQHE of ®CEs in a negative B* becomes essen¬ 
tial when we consider the spin degree of freedom, because 
it is the only known way to explain the non-fully spin- 
polarized EQHE states here, e.g. the spin singlet state 
at z^ = 2/3. (Recall that for spinful states, p-h symmetry 
relates z/ t o 2 — u.) An extensive experimentaP^^^ and 
theoreticaP^^^^^^^®^ literature on spin phase transitions 
has validated the explanation of the z/ = (n + 1)/(2n + 1) 
as z/* = n + 1 IQHE of ^CEs. 

The validity of '^ri/( 2 n±i) foi* fhe incompressible 
states has been established by extensive numerical 
studiej^ ^^ l ^^ l ^^^ tnni, We will make the assumption that 
^n/(2n±i) remain valid to arbitrarily high n, i.e., that 
the compressible region around z/ = 1/2 consists of un¬ 
resolved IQHE states of composite fermions. We stress 
that we cannot rule out the possibility that the ®CEES 
and ^CEES are in reality topologically distinct, as pro¬ 
posed in Ref. [26l and a spontaneous breaking of the p-h 
symmetry selects one of them. This would happen, for 
example, if the half filled ground state were unstable to 
CE pairin^^IH^ and ®CEES and ^CEES are the normal 
states of the topologically distinct Moore-Read Pfafiian 
and anti-Pfafhan paired-CE state^^^^^. Nonetheless, the 
presently known facts do admit the possibility of a p- 
h symmetric CEES, and our aim here is to deduce its 
properties, so experiments may distinguish between the 
different proposals. 

We have calculated the pair correlation function for 
V = nl(2n + 1) up to 7/15 using the Metropolis Monte 
Carlo method. Eor technical reasons, we find it con¬ 
venient to use the standard spherical geometr}!^ see 
SM for detailJ^. The results extrapolated to the ther¬ 
modynamic limit apply to the planar geometry of the 
experiments. The spherical analogs of the above wave 



FIG. 1: (a) Pair correlation function g{r) as a function of r/^, 
where r is the arc distance on the sphere. The projected wave 
functions in Eq.j^have been used for its evaluation. The solid 
lines are fits using Eq. [^for the initial oscillations. For clarity, 
the curves (except for 5/11) have been shifted up or down by 
multiples of 0.02. (b) The calculated thermodynamic values 
of k^i as a function oi n. The mean-field values and 
vW — n) are also shown for reference. 


functions, as well as the details o f LLL projection can 
be found in the hteratur (j^^FQFi [ wave functions 

considered below have uniform density and are trans- 
lationally invariant (i.e. have orbital angular momen¬ 
tum L = 0 on the sphere). The g{rys for the largest 
systems in our study are shown in Fig. [^a). For incom¬ 
pressible systems the pair correlation function is expected 
to decay in a gaussian manner in the limit of r ^ oo, 
but there is a range of intermediate r where it exhibits 
well defined oscillations from which a Fermi wave vec¬ 
tor can be extracted. In fitting g{r) to Eq. we avoid 
very small r (where short distance correlations are im¬ 
portant) and very large r (where curvature effects be¬ 
come non-negligible). Erom the results for finite systems 
we obtain the thermodynamic limits for the k^i (se^^. 
We find that very large systems (A^ > 100) are needed 
for a satisfactory thermodynamic extrapolation of k^i. 
The thermodynamic limits are shown in Eig. [^b). [We 
have assumed exact p-h symmetry, which implies that 
the k^i at z^ = (n + l)/(2n + 1) is the same as that at 
u = njiyin + 1).] The range of k^£ includes uncertainly 
in the fits (estimated by linear and quadratic fitting in 
1/A^ for g{r)) as well as uncertainty due to the curvature 
of the spherical geometry (estimated by considering fits 
with r chosen as the chord or the arc distance). Eor ref¬ 
erence, Eig. yjb) also shows the values k^^^i = ^/^ and 
_ ^2(1 — z/) as expected from a ME picture for 
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FIG. 2: Same as in Fig. but for the unprojected wave 
functions 4^n/(2n±i)- Also shown for reference is = 

\/^ corresponding to the Chern-Simons mean held theory 
for ^CFs. 

the ^CFFS and ^CFFS. 

For u = 1/2, we have estimated k^i by an extrapola¬ 
tion, in the spherical geometry, of hlled shell CF systems 
at zero effective flux^ occurring at N = v?. For techni¬ 
cal reasons, we are not able to go to systems with > 81 
(which requires hlling the 10th Landau-like level of com¬ 
posite fermions, where the numerics become unstable). 
We have therefo re also studied the CFFS in the torus 
geometr} ^^^l^^ l ^^ l where we can go up to = 153, and 
hnd that the results are consistent with our spherical re¬ 
sults. Results in the torus geometry are presented irP^. 

For V away from 1/2, our calculated k^^ is close, but 
not equal, to the smaller of and ^^2(1 — v). Both 
from extrapolation of the results from the sequences 
nj(2n ± 1) and from calculations directly at v = 1/2, 
our calculations suggest that the CFFS area at = 1/2 
slightly deviates from the value expected from the Lut- 
tinger rule. 

The physics of the CFFS at v = 3/2 is analogous to 
that at z/ = 1/2 once the B dependence of the density of 
either® CFs or ^CFs in the spin-reversed LL is accounted 
fQ j23 | 74 | jy — 2 / 4 ^ both n/(4n+l) and n/(4n —1) are 

understood only in terms of ®CFs, and thus one expects 
k^i ^ ^/^ (as opp osed t o kpi ^ ^^2(1 — z^)), as ob¬ 
served experimentall}^^^^ and also in our calculation^^. 
Analogous consideration for the CFFS at u = 3/4 gives 
k*p£ « V2(l-i/)- 

We next investigate how robust the CFFS area is to 
LL mixing. LL mixing requires a formulation in terms 
of electrons (rather than holes of the LLL), and the 
LLL electronic wave functions in Eq. can be used as 


a starting point to address this issu^^. A realistic treat¬ 
ment of LL mixing is outside the scope of the current 
study, but we have considered the “unprojected” Jain 

wave functions '^Z/(2n±i) = which 

contain some amplitude outside of the Even 

though they do not give a realistic account of LL mix¬ 
ing, it is likely that they are adiabatically connected 
to the projected wave functions (as explicitly demon¬ 
strated for V = 2/d^, and hence to the actual Coulomb 
ground states. For these wave functions, the 5 '(r)’s of 
V = nl{2n — 1) and v = nji^n + 1) are identical for a 
given N when plotted in units of the sphere radius, which 
in the thermodynamic limit implies the relation 

/2n + 1 \ 

= ( 9 ) 

The calculated values of k^'^^ using the ®CF description 
(see SM for details) are shown in Fig. Our calcula¬ 
tions thus provide evidence that k^^^ depends on LL mix¬ 
ing. Another approximate wave function with LL mixing 
is the CS ME state “ Zk)/\zj - 

^/c|]^^±n(^*)- Given that its g{r) coincides with that for 
^±n(^*): we get i.e., k^^I = \/^ for 

all V = nj(2n ± 1). For the unprojected or the CS-MF 
wave functions, /cp^ does not obey particle-hole symme¬ 
try, as expected in the presence of LL mixing. 

In summary, we have shown that, within the micro¬ 
scopic theory of composite fermions, it is valid to consider 
electron (or hole) based composite fermions for v < 1/2 
as well as v > 1/2. We have calculated the CF Fermi 
wave vector in the vicinity of v = 1/2 and find that it 
is closer, but not equal, to the smaller of y/47rpe and 
y/47rph- la terms of electron based composite fermions, 
this implies that the Luttinger theorem is slightly (sub¬ 
stantially) violated for z^ < 1/2 (z^ > 1/2). At z^ = 1/2, 
our results suggest, but do not conclusively demonstrate, 
that the k^ differs slightly (by a few percent) from the 
value predicted by the Luttinger area rule. We 

also provide evidence that k^ varies as a function of LL 
mixing. 

Note Added: Since the completion of this work, several 
other articles have appeared addres sing th e nature of the 
CFFS and the role of p-h symmetr}ffSj^. 
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Supplemental material 

The Supplemental Material (SM) sections contains (i) 
a primer on the spherical geometry; (ii) a discussion 
of particle-hole (p-h) symmetry within the composite 
fermion (CF) theory; (iii) the relation between pair cor¬ 
relation functions of two states related by the p-h sym¬ 
metry for finite N; (iv) extrapolations of the CF Fermi 
wave vector and the power law exponent a; (v) a discus¬ 
sion of the CF Fermi Sea (CFFS) at u = 1/4; and (vi) 
our calculations for the CFFS at = 1/2 and z/ = 1/4 in 
the torus geometry. 

A note on the notation: Quantities defined with the 
subscript or superscript ‘e’ are for electrons and those 
with superscript or subscript ‘h’ are for holes. 


I. HALDANE’S SPHERICAL GEOMETRY 


For many of our calculations we employ Haldane’s 
spherical geometr}!^ shown in Fig. SI In this geom¬ 
etry, N electrons move on the surface of the sphere and 
the magnetic field perpendicular to the surface is gener¬ 
ated by a monopole of strength Q placed at the center 
of the sphere. The total flux through the surface of the 
sphere is 2Qhc/e and the radius of the sphere is ^/Qi. 
The electron coordinates are denoted by ft = {O^cj))^ and 
the single particle eigenstates, called monopole harmon¬ 
ics, are given by (see, e.g. the review in Ref. [12]): 




xE 


l—m 

s=0 




(SI) 

{u^uY 


{21 + !)(/ — m)\{l + m)\ 
" ^ [ 4^ (/-Q)!(/ + Q)! 

u = cos( 6 >/ 2 )e *2 

V = sin(^/2)e“*2 


(5 2 ) 

(53) 

(54) 


where the orbital angular momentum I takes values 
IQ/IQI + 1, |Q| + 2 ,... corresponding to the n = 0 
(lowest), 1 (second), 2 (third),... Landau Levels (LLs) 
and m = ——I + 1 ,..., / — 1 , /. Consequently the degen¬ 
eracy of each LL is 2/ + 1. Positive and negative values 
of Q respectively stand for a flux extending radially out¬ 
ward and inward through the surface of the sphere. The 
wave functions in the lowest Landau Level are obtained 
by setting / = Q in the above equation which gives: 


^Q,Q,m 


- 2 Q-M f 2Q \ 
. 4:71 \ 2(3 — m) 




(S5) 

The wave function for a state with an integer number of 
completely filled LLs is a Slater determinant formed from 
the filled single particle states. Using the above equation, 
one can show that the form of the wave function for one 


filled LL (apart from a normalization) has the following 
nice form: 


^1 = - UkVj) (S 6 ) 

j<k 

In the main article, the wave functions of composite 
fermions have been written in coordinates appropriate 
for the planar geometry. They can be expressed in the 
spherical geometry in the standard fashion explained in 
the literatur^i^. We have also used the torus geometry 
(see below) for an additional substantiation of some of 
our results. The choice of which geometry is used is a 
matter of technical convenience. For example, the spher¬ 
ical geometry is more suitable for the FQHE states at 
n/(2n ± 1) than the other two geometries. We carefully 
extrapolate all of our results to the thermodynamic limit 
to eliminate any effects of the geometry employed. 



FIG. SI: A schematic of the Haldane sphere, in which N 
electrons reside on the surface of sphere. A Dirac monopole 
at the center produces a flux of 2Q (an integer) in units of 
the flux quantum 0o = hc/e. 


H. PARTICLE HOLE SYMMETRY IN THE CF 
THEORY 

In the CF theory, the states at z^ = n/(2n + 1) maps 
into z/* = n of ®CFs, whereas its hole partner at z^ = 
1 — nf{2n + 1 ) = (n + 1 )/ {2n + 1 ) maps into z/* = n + 1 
of ^CFs. This may seem to be inconsistent with p-h 
symmetry, but we show that is not the case. 

From the standard CF theory on sphere, the state at 
V = nj{2n + 1 ) occurs at flux value 

-I- 1 

2Q = - (2 + n) (S7) 

n 

Using the relation 2(5* = 2Q — 2Nq -h 2 we have 2(5* = 
NQ/n — n, which corresponds to precisely n filled A levels 
of ®CFs (CF Landau-like levels are termed A levels). The 
single particle orbital angular momentum of the topmost 
occupied AL is 


/: = Q:+n-i 


(S 8 ) 














Now let us ask what happens if we model this state in 
terms of ^CFs. In terms of A^h = 2(5 + 1 — A^e holes, the 
V = (n + 1)/ (2n + 1) state occurs at 

2t7 -4- 1 

2Q = - (1 - «) (S9) 

n + 1 

This maps into a negative effective flux of 2(5^ = 
—N\^/{n + 1) + n + 1 for ^CFs, which corresponds to 
n + 1 filled ALs of ^CFs. The single particle orbital an¬ 
gular momentum of the topmost occupied AL is 

ll = \Ql\+n. (SIO) 

Using the above relations, it follows that /* = 
topmost fully occupied AL has the same single particle 
orbital angular momentum, and thus the same number of 
particles, independent of whether we model the problem 
in terms of ®CFs or ^CFs. This implies that the ®CFs 
and ^CFs descriptions produce identical structures, at 
the mean field level, for the excited bands for energies 
up to nhuj^^ where huj^ is the effective cyclotron energy 
of composite fermions. We further note that the cy¬ 
clotron energies for both are given by fkj* = fie5*/m*c = 
fie5/(2n + l)m*c, where B* is the effective magnetic field 
for composite fermions and m* is their mass. This im¬ 
plies that that for a given magnetic field, the two de¬ 
scriptions produce the same cyclotron energy, and thus 
the same excitation gaps at the mean field level. As noted 
in the text, the microscopic wave functions of composite 
fermions also obey the particle-hole (p-h) symmetry to a 
very good approximation. 


III. RELATION BETWEEN AND IN THE 
SPHERICAL GEOMETRY 

In this section we shall derive the exact formula re¬ 
lating the pair correlation functions of two states re¬ 
lated by particle-hole symmetry, in the spherical geom¬ 
etry for finite N. (The reiation derived in the main 
text assumes disk geometry and transiationai invariance, 
and thus is vaiid oniy in the thermodynamic limit.) 
We then test how well this relation is satisfied by the 


pair correiation functions of T^lll Ylj<k(^j ~ and 

Vlll Iljcki^j ~ • That wouid be a test of the 

degree to which these wave functions are reiated by p-h 
conjugation. 

The iowest Landau ievei projected eiectron fieid an¬ 
nihilation operator in the spherical geometry is written 
as: 


m=-\-Q 

4'(r)= ^ c^YQ,Q,m(r} (Sll) 

m= — Q 

where r = fl = {u^v) = (6>, 0) is the coordinate on the 
sphere [0 and 0 are the polar and azimuthal angles on the 
sphere, u and v are spinorial coodinates defined as: u = 
cos(6>/2)e*'^/^ and v = sin(6>/2)e“*'^/^] and are the 

so-calied sphericai monopoie harmonic^^. The foiiowing 
identities are satisfied by the eiectron fieid operator: 


{T^(r), T(r')} = S{r',r) 
{(l)^\^{r)^\r')\(t)^) = {1- X)6{r,r'), 

= A5(r',r) (S12) 

wher^^ 


(5(r, r') 


2Q H- 1 
47rQ£‘^ 




N 

2Q + 1 


(S13) 


To derive the reiation between the pair correiation 
functions at u and 1 — u we start with Eq. 5 of the 
main text. Using the reiations in Eq. S12 we can show 
that: 


gi-u{r,r') = - 6{r,r'f] 

Pl-y 

+ ^g,{r,r') (S14) 

Pl-V 

Next we make use of the fact that incompressibie states 
have both rotationai and transiationai invariance, and 
choose r' to be the north poie. Substituting the vaiues 
from Eq. |S13[ we find that 


9i-v{rc) = 


(2g + l)(2Q + l-2Ar) 
(2Q + 1 - NY 


1 - 


1 - 


1 2Qx 

(5^)1 ) + ( 20+4-jv )'^-<’'°> 


2Q + 


(S15) 


where Vc is the chord distance of the point r from the 
north poie. It can be verified that this reduces to Eq. 7 
of the main text in the thermodynamic limit A/", 2Q oc 
with N/2Q u. 

In Eig. I^we show a comparison of the pair correla¬ 


tion functions of ^(n+i)/(2n+i) and (^p-hT^/(2n+i)- (The 
g{r) of the lat ter is obtained from the g{r) of T^/( 2 n+i) 


following Eq. S15). The excellent agreement between 


the two pair correlation functions confirms that the wave 
functions '^ri/(2n-\-i) and are reiated by p-h 
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symmetry to an extremely good approximation. 


IV. EXTRAPOLATIONS IN THE SPHERICAL 
GEOMETRY 


In this section, we provide some further details associ¬ 
ated with the results shown in the main article. Fig. 
(S4) shows the extrapolations of to the ther¬ 

modynamic limit for the projected (unprojected) wave 
functions of Eq. 8 of the main article. The extrapolated 
values are shown in Fig. 1 and Fig. 2 of the main arti¬ 
cle. In each case, we ascertain the range of by fitting 
it as a function of 1/V with a linear and a quadratic 
form, using both the arc and the chord distances; we be¬ 
lieve the error estimated in this fashion adequately takes 
into account finite size as well as curvature effects. We 
also show the extrapolation of the power law exponent 
a as defined in Eq. 1 of the main text using the 

projected (unprojected) Jain wave functions at various 
filling factors in Fig. [S^(S6). 


V. CF FERMI SEA AT = 1/4 

The FQHE states near = 1/2 can be obtained in 
two seemingly distinct fashions. For example, the state 
at n/(2n — 1) can be thought of as either filling factor 
z/* = n of ®CFs in a negative effective magnetic field, or 
as the electron partner of z/* = n — 1 of ^CFs in a positive 
effective magnetic field. For the FQHE states near u = 
1/4 the situation is different. The only way to obtain 
states at z^ = nji^n + 1) (z^ = nj{^n — 1)) is in terms of 
®CFs at z/* = n in positive (negative) effective magnetic 
field. As a result, we have a unique CFFS at v = 1/4 
and the apparent dichotomy encountered above near v = 
1/2 does not present itself near v = 1/4. Viewing this 
state in terms of holes at z/h = 3/4, and then composite 
fermionizing the holes does not produce a CFFS of ^CFs 
at zero effective magnetic field. From this perspective, it 
is to be expected that the the Fermi wave vector will be 
close to /cp-^ = y47rpe near v = 1/4. 

We have performed extensive calculations to confirm 
this expectation. Our calculations are carried out at fill¬ 
ing factors V = n/(4n±l) which arise from IQHE of ®CFs 
carrying four vortices. The projected Jain wave functions 
for the ground states at these filling factors are given by: 


For V = 1/4, akin to the case of z^ = 1/2, we have calcu¬ 
lated /cp by considering filled shell CF systems occurring 
at zero effective flux. 

Our results in Fig. S7 (S8) show the extrapolations 
of /cp^ {k^^t) to the thermodynamic limit for the pro¬ 
jected (unprojected) wave functions. The extrapolated 
values are shown in the bottom-most panels of the re¬ 
spective figures. To obtain these plots we follow exactly 
the same procedure as illustrated above for the sequence 
n/ (2n±l). We also show the extrapolations for the power 
law exponent a as defined in Eq. 1 of the main 

text, obtained from the projected (unprojected) wave 
functions at various filling factor s of c omposite fermions 
carrying four vortices in Fig. S9 (S1Q|). 


We find that the obtained Fermi wave vector in the 


vicinity of z/ = 1/4 is close to the value for both 

the projected and unprojected wave functions. At pre¬ 
cisely z^ = 1/4 we find that the Fermi wave vector ob¬ 
tained from the projected wave functions is close to but 
slightly different from the mean field value of 1/\/2^, sug¬ 
gesting a slight violation of the Luttinger theorem. For 
the unprojected state at 1/4 we find that the Fermi wave 
vector is consistent with the mean field value of 

The situation close to u = 3/4 is analogous. The 
CF Fermi wave vector will be close to k^£ = y47rph = 


Y^2(1 — u) for situations where p-h symmetry is valid. 

Since we are not able to go beyond 81 particles (81 par¬ 
ticles) for u = 1/2 (z/ = 1/4) in the spherical geometry, 
we have also considered the CFFS in the torus geometry, 
where we have studied the CFFS for up to 153 particles 
(121 particles). The calculation in the torus geometry is 
described next. 


VI. CF FERMI WAVE VECTOR IN THE TORUS 
GEOMETRY 

A brief outline for the calculation on the torus geom¬ 
etry for u = l/2p, where p is an integer, is as follows: 
We consider a periodic geometry in the plane spanned 
by two nonparallel vectors, Li an L 2 . Using the sym¬ 
metric gauge A = (—H^/2, Hx/2,0) for a homogeneous 
magnetic field B = (0,0,H), the wave functions in the 
lowest Landau level have the form 

tp{x,y) = exp f{z), (S18) 


N 

^ra/(4ra±l) = "PlLL H {Zj “ ZkY^±n (S16) 

j<k=l 

where Plll denotes LLL projection and is the wave 
function of n filled LLs, with = [4>^]*. The un¬ 
projected Jain wave functions for the ground states are 
defined as: 

N 

^n/(4n±l) {^j ^±n 

j<k=l 


where z = {x iy)l\f2£ and f{z) is an entire function. 
Consequently, the many-body states have the form 

= exp F{zi,Z 2 ,...,zn), (S19) 

where F is entire. We will refer to the first factor as the 
nonanalytic part. We require periodic boundary condi¬ 
tions in each coordinate,^ 

= exp{i(pj)'^, 


(S17) 


(S20) 
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P8 01i S 0( W f)0 8 il0MI0 1 l^^ 


® breverse flux: N = 75, z/ = 3/5 CFj 
•from p-h of: N = 52, z/ = 2/5 CF| 

'0 2 4 6 8 10 12 14 

r/£, r is the chord distance 


bC 



4 / 

9 



oreverse flux: N = 60, f = 5/9 CF 
• from p-h of: N = 52, v - 4/9 CF 


r/£, r is the chord distance 




FIG. S2: Pair correlation function g(r) as a function of r/^, where r is the chord distance. The blue circles give the g(r) 
obtained directly from the lowest Landau level projected wave functions T(^+i)/( 2 n+i) foi* 3/5 (top left panel), 4/7 (top right 
panel), 5/9 (bottom left panel) and 6/11 (bottom right panel). The red dots give the g(r) for Cp-h4^n/(2n+i)? namely the hole 
conjugates of projected wave functions T^/( 2 n+i)- 


where j G {1,2}, (j)j are fixed constants, and t^(L) is 
the magnetic translation operator acting on the com¬ 
plex coordinate defined in terms of the Cartesian ones 


t^(L) = exp 


L • (V - {ie/h)A) - —{Lx Tm)' 


The two conditions in Eq. (S20) are compatible only if 

|Li xL2\ = 2Tre^N^, 


(521) 
lyif 

(522) 


where is an integer; i.e., the parallelogram spanned 
by Li and L 2 encloses an integer number of flux quanta. 

The unprojected CF Fermi sea is simply the product 
of a Slater determinant of plane waves and the Jastrow 
factor that binds an even number 2p of flux quanta of the 
many-body wave function to each electron: 

Det(e**'^--^)vl>i/2p, (S23) 

where '^i/ 2 p is the bosonic Laughlin at filling 

factor TJ = l/2p, which is a function oi N = N(f)/2p elec¬ 
tron coordinates, and {kj} are the momenta of the filled 
CF states, taken from reciprocal lattice that corresponds 
to the direct lattice spanned by Li and L 2 . 

When the lowest Landau level projection is imple¬ 
mented in the periodic geometry, differentiation with re¬ 
spect to Zm becomes a (magnetic) single-body translation 


operation. In particular, the CF Fermi sea state i^ 

Det(ti(-d^))^i/2p, (S24) 

Here, d = (mLi + nL 2 )/A^(/), and the momenta of the 
filled CF states are generated as k = {dyj 
With t he d eterminant spelled out, there are N\ terms 
in Eq. (S24), making it computationally intractable. We 


will use the convenient approximate projection scheme 
elaborated by Shao et This leads to the wave func¬ 
tion 


^CFFS = Detij 


e^>^l[a^P{zi-Zk-2{dj-d)) 




]{zi - d) exp 



(S25) 


where dj = {dj^^ + idj^y)fV2i, d = {1/N) Y.f=i dj, -Fc.m. 
in an analytic center-of-mass wave function, and (7{z) is 
the Weierstass’ sigma function appropriate to the lat¬ 
tice spanned by the complex numbers Li = {Li^^ + 
iLx^y)l\f2i and = {L^^x + iL 2 ,y)lV 2 t 


a{z) 


di{'Kz/Li\T) / i{'Kz/Li)'^ \ 

(7r/Li)i?;(0|T) V7r('r - 'T*) ) 


(S26) 


with T = L 2 /L 1 . Using the properties of the a function, 
one obtains by elementary steps that the periodic bound- 
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FIG. S3: Thermodynamic extrapolation of the Fermi wave vector for the projected Jain wave function at various Filing 
factors. In this, as well as the following seven figures, the empty (filled) symbols correspond to the values obtained from the 
chord (arc) distance on the sphere and the thick (dashed) lines show linear (quadratic) hts to these values as a function of 1/N. 


ary conditions in Eq. (S20) imply the following require¬ 
ments {j G {1,2}) for the center-of-mass wave function: 


Fc.mX^ + ^ 2 ) 

Fc.mXZ) 


= exp ( i(l)j + 


p\Lj 




+ 


2pL* 


(S27) 


As the pair correlation should not depend on the center- 
of-mass motion, we are content with finding a particular 
solution to Eq. (S27). Seeking it in the form that is 
structurally analogous to the a function, 

f7rf3Z 


Fc.m. {Z) = exp 


faZ^\ ^ fTr(3Z 

lTrj"'(Tr 


|/3r 


(S28) 


the boundary conditions fix /3 = 2p and a = if 

01 = 0 and 02 = tt. 


Specializing now for a square unit cell, Li • L 2 = 0 
and |Li| = IL 2 I, the pair correlation was obtained up to 
N < 153 at u = 1/2 and up to < 121 at u = 1/4. 
The functional form in Eq. (1) of the paper was fitted 
to obtain the Eermi wave vector. Examples are shown in 
Eig. [Sill The e xtrap olation to the thermodynamic limit 
is shown in Eig. S12 At = 1/2, we get k^l = 1.03(1) 
for a linear fit in 1/N. When we allow both linear and 
quadratic fits, the range is given by k^l = 1.02 — 1.07, 
which is consistent with the results from the spherical 
geometry quoted in the main text. At u = 1/4, we 
get k^l = 0.66(2) for a linear fit; with both linear and 
quadratic fits we get the range k^l = 0.61 — 0.69, which 
is again consistent with the results from the spherical 
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FIG. S4: Thermodynamic extrapolation of the Fermi wave vector for the unprojected wave functions 4^n/(2n±i) defined 

in the main text at various filling factors of the form h' — n/(2n+ 1). 


geometry. 

A note regarding the scatter in our data is in order. 
The main problem of CF Fermi sea on the torus is that 
one has to select the wave vectors from the reciprocal 
space (i.e., wave vectors compatible with the boundary 
condition) that are occupied. Shao et propose that 
the choice that gives the lowest energy Fermi sea for non¬ 
interacting particles, which is a roughly circular disk in 
k-space, is appropriate. However, for relatively smaller 
system sizes, the CF Fermi sea is not quite circular, which 


is a source of finite size effects. 

We have performed a similar analysis for the unpro¬ 
jected CF Fermi sea. The results are shown in Fig. |S14| 
for v = 1/2. For a linear fit we find k^l = 0.96(1). 
When we allow both linear and quadratic fits, the range 
is given by k^l = 0.96 — 1.00, which again is consistent 
with the range quoted in the main text. At = 1/4, we 
get k^l = 0.678(2) by linear extrapolation, or the range 
k^l = 0.67 — 0.69 with both linear and quadratic fits. 
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FIG. S5: Thermodynamic extrapolation of the power law exponent a (defined in Eq. 1 of the main text) evaluated using the 
projected Jain wave functions at n/(2n ± 1). The bottom-most panel shows thermodynamic values of a as a function of the 
Filing factor ly. 
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FIG. S6: Thermodynamic extrapolations for the power law exponent (as dehned in Eq. 1 of the main text) evaluated 
using the unprojected wave functions 4^n/(2n+i) ^^r Filing factors of the form nl{2n + 1). The bottom-most panel shows the 
thermodynamic value of as a function of the Filing factor ly. 
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FIG. S7: Thermodynamic extrapolation of the Fermi wave vector k^i obtained using the projected wave functions '^n/{ 4 n±i) 
at nl{An ± 1). The bottom-right most panel shows thermodynamic values of kpi as a function of the filling factor ix. 



























































































16 



0.601(2) 

0.6 

0.59(1) p 0.59(1) 



0.75 “ O Chord distance 

— Linear, X = 1-6 
-- Quadratic, X =1.6 
□ Arc distance 
0.7 h Linear, = 1-5 
- Quadratic, x^^ed = ^ 

0.65 


0.595(3) 


V = 3/13 

-ffl ■■ ■'™~ ■' 'ffi 


0.8 


O Chord distance 
— Linear, X ^ed = ^ 

- - Quadratic, x^^ed = ^ 

□ Arc distance 
Linear, X ^ed = 

0.75 h “ Quadratic, x^red = 0-9 




0.7 


0.662(1) 

0.65 7 
0.645(7) 



0.01 

O Chord distance 
— Linear, X ^ed = 

-- Quadratic, x^^ed = 
□ Arc distance 

Linear, x^red = 

- Quadratic, x^red = 


0.02 

1/N 


0.03 



V = 4/17 


0.02 

1/N 


- 0.604(7) 


0.01 


0.02 

1/N 


0.03 


O Chord distance 




O Chord distance 

— Linear, x^red = 2-2 


0.8 


— Linear, x^red = ^7 

- - Quadratic, X red “ ^ 



-- Quadratic, x^red = 

□ Arc distance 




□ Arc distance 

Linear, x a = 2.5 

’ red 




Linear, x red = 

■ - ■ Quadratic, x red “ 


-- 0.75 


- ■ Quadratic, X^red “ 



PLh 

4 ^ 


0.65 

0.635(3) , 
0.63lU) 
0.625(4) 



0.005 


0.01 0.015 

1/N 




FIG. S8: Thermodynamic extrapolation of the Fermi wave vector obtained using the unprojected wave functions ^^/(dnii) • 
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FIG. S9: Thermodynamic extrapolation of the power law exponent a obtained using the projected Jain wave functions at 
Filings nl{4:n ± 1). The bottom-right most panel shows thermodynamic values of ct as a function of the filling factor jy. 
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FIG. SIO: Thermodynamic extrapolation of the power law exponent obtained using the unprojected wave functions 
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FIG. Sll: Fit of the Friedel oscillations in the pair correlation function of the composite fermion Fermi sea state at h' — 1/2 
(left) and v — \/A (right) in the torus geometry. The results shown are for the largest systems considered in this work. 




FIG. S12: The Fermi wave vector (left panel) and the power law exponent a (right panel) of the projected composite fermion 
Fermi sea state at z/ = 1/2 on the torus geometry as a function of 1/N. The linear extrapolation to the thermodynamic limit 
is also shown. 
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FIG. S13: The Fermi wave vector (left panel) and the power law exponent a (right panel) of the projected composite fermion 
Fermi sea state at z/ = 1/4 on the torus geometry as a function of 1/N. We show both the linear and the quadratic extrapolations 
to the thermodynamic limit. 




FIG. S14: The Fermi wave vector (left panel) and the power law exponent (right panel) of the unprojected composite 
fermion Fermi sea state at z/ = 1/2 on the torus geometry as a function of 1/iV. The linear extrapolation to the thermodynamic 
limit is also shown. 
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FIG. SI5: The Fermi wave vector (left panel) and the power law exponent (right panel) of the unprojected composite 
fermion Fermi sea state at z/ = 1/4 on the torus geometry as a function of 1/N. We show both the linear and the quadratic 
extrapolations to the thermodynamic limit. 
















































